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Effective field theories of two-dimensional lattice models of fluctuating loops are constructed 
by mapping them onto random surfaces whose large scale fluctuations are described by a Liouville 
field theory. This provides a geometrical view of conformal invariance in two-dimensional critical 
phenomena and a method for calculating critical properties of loop models exactly. As an application 
of the method, the conformal charge and critical exponents for two mutually excluding Hamiltonian 
walks on the square lattice are calculated. 
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The study of critical phenomena in two dimensions has 
led to a rather remarkable interplay between statistical 
mechanics and field theory. The picture that has emerged 
portrays critical fluctuations in a variety of systems as be- 
ing described by a conformal field theory (CFT), i.e., a 
quantum field theory invariant under local scale (confor- 
mal) transformations. This view has proven very fruitful 
as it has answered many questions pertaining to two- 
dimensional phase transitions. For example, by study- 
ing the representations of the Virasoro algebra, which is 
the algebra of conformal transformations, we have come 
to understand why critical exponents in two-dimensions 
are typically rational numbers jlj . Lists of possible crit- 
ical exponents appear here in very much the same way 
one discovers the quantization of angular momentum by 
studying the representations of the algebra of rotations. 

Progress in understanding two-dimensional critical 
phenomena using CFT's has largely come about by clas- 
sifying these field theories and associating them with the 
scaling limits of specific microscopic lattice models 
Still a fundamental practical question remains: Given a 
lattice model how does one go about constructing a con- 
formal field theory of its scaling limit? It is this question, 
limited in scope to lattice models of fluctuating loops 
H, which is addressed here. Since most canonical two- 
dimensional models can be recast as loop models [B , the 
results we find are general in nature, and provide a ge- 
ometrical view of criticality and conformal invariance in 
two-dimensions. In particular, we find that CFT's of 
loop models can be constructed explicitly via a mapping 
of the loop model to a model of a fluctuating surface. 
The basic idea is to think of loops as contour lines of 
a random surface; the field theory which describes the 
large-scale fluctuations of this surface is a conformally 
invariant Liouville field theory. 

Liouville field theory provides a general framework for 
calculating correlation functions in CFT's; this is the well 
known Coulomb gas representation pioneered by Dot- 
senko and Fateev ||. The central result of this paper 
is the discovery that different terms in the Liouville ac- 
tion, some of which were originally introduced on for- 
mal grounds, have a concrete geometrical interpretation 
in loop models. Furthermore, the coupling constants, 
most notably the effective stiffness of the surface K, are 



fixed by geometry, thus providing a method for calculat- 
ing critical properties of loop models exactly. To illus- 
trate this we calculate the conformal charge and critical 
exponents for a new loop model defined on the square 
lattice. For vanishing loop weights this model defines a 
new universality class of Hamiltonian walks; these are 
self-avoiding random walks that visit all the sites of the 
lattice. Hamiltonian walks have been used in the past to 
model configurational statistics of polymer melts [|| . 

a. Loop models Given a two-dimensional lattice C a 
loop configuration Q is defined by a set of loops, which are 
closed self-avoiding walks along the bonds of L. The par- 
tition function is Z = ^Zg ex P[~/?-ff (<?)], where (3H(Q) is 
the reduced loop-Hamiltonian. 

As remarked in the introduction, many two- 
dimensional models can be recast as loop models. For 
example, in the Potts model one arrives at a loop model 
description by considering the graphs generated by the 
high-temperature expansion of the partition function [|J ; 
in the O(n) model, the partition function takes the form 
of a loop model once integration over the spin degrees of 
freedom is carried out fflj. Here we focus our attention 
on fully packed loop (FPL) models whose allowed config- 
urations Q satisfy the constraint that every vertex of the 
lattice is visited by a loop, and 

PH{G) = Y^^ N M- (1) 

is the loop Hamiltonian. Here we have allowed for rif 
flavors of loops; = exp(— /ij is the loop weight, and 
Ni is the number of loops of flavor i. In the allowed con- 
figurations Q only different flavored loops can cross. An 
important example of an FPL model is provided by the 
critical Q-state Potts model, which can be mapped to an 
FPL model on the (oriented) square lattice with n r — 1 
and n = y/Q 0. Recent interest in FPL models Q [ll]] 
has been sparked by their close relation to the problem 
of Hamiltonian walks: from the Bethe ansatz solution of 
the FPL model on the honeycomb lattice exact critical 
properties of Hamiltonian walks on a non-oriented lattice 
were obtained for the first time ||. 

In order to illustrate the general relation between loop 
models and Liouville field theories, we calculate the exact 
conformal charge and critical exponents for the two-flavor 
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(rif = 2) fully packed loop model (FPL 2 ) on the square 
lattice, with loop weights rt\ — n 2 = n > 0. The two loop 
flavors of the model will be denoted black and grey, in 
accord with Fig. ^a; the possible vertex configurations of 
the FPL 2 model are shown in Fig. |l|b. The allowed loop 
configurations are the same as in the dimer-loop model 
JL^I ) and the FPL model on the square lattice jy]], the 
difference being in the loop weights. For the dimer loop 
model n\ = 2, 112 = 1, while the FPL model generalizes 
this to m = n > 0, n% = 1. Obviously the FPL and 
FPL 2 models share a common point: n = 1. We will 
show later that for n = 1 there is excellent agreement 
between the numerical results of Rcf. [O and our exact 



expressions. 
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FIG. 1. The two-flavor fully packed loop (FPL 2 ) model on 
the square lattice: a) loop configuration Q; b) allowed vertex 
configurations; c) two vertex configurations of oriented loops 
(of 24 total) and their weights - s = exp(i7reo/4); d) mapping 
to the coloring model - "e" and "o" denote sites on the even 
and odd sublattice of the square lattice. 

It is instructive to look at the n — ► (/i — * 00) and 
n — > 00 (fi — > —00) limits of the FPL 2 model. From 
Eq. (|l]) we infer that in the first limit configurations with 
one black and one gray loop are selected. These self- 
avoiding walks are Hamiltonian and they can be thought 
of as a model of two mutually repelling species of poly- 
mers in the melt. Later we will show that this model 
defines a new universality class of Hamiltonian walks. In 
the limit n — > 00 the selected loop configurations maxi- 
mize the number of loops, each being of minimum length 
four. In this case the model has a finite correlation length 
(£), measured by the average loop size, while in the case 
of Hamiltonian walks the correlation length obviously di- 
verges with system size. Later we will argue that at n = 2 
a transition occurs between a critical phase (£ — ► 00) and 
a long range ordered one (finite £). 

b. Liouville field theory To calculate the critical prop- 
erties of a loop model we construct a random surface for 
which the loops are contour lines. This construction can 
be broken up into three steps, which are illustrated here 
for the FPL 2 model. 

First we orient all the loops; for every loop configu- 



ration Q this gives 

2 JV 1 (<?)+jV 2 (S) oriented-loop confi gura- 
tions. In order to recover the proper loop weights (n) 
every left /right turn of the oriented black and grey loops 
is assigned a weight A = exp(±i7reo/4); if the loop does 
not make a turn, A = I is the appropriate weight (see 
Fig. |l]c). On the square lattice the difference between the 
number of left and right turns that a closed loop makes is 
±4, so the A's along a loop will combine to give a weight 
of exp(±i7reo), depending on the loop orientation. After 
summing over the two possible orientations, the correct 
loop weight 



2cos(7reo) 



(2) 



is recovered. 

There is an important caveat to the above construc- 
tion. Namely, if we define the loop model on a cylin- 
der, then the loops that wind around the cylinder are 
weighted incorrectly, with n — 2, since the number of 
left and right turns for these loops is equal. To correct 
this a boundary term is introduced in the effective action 
for the FPL 2 model - S b in Eq. (g) §. 

Second, we map configurations of oriented loops to con- 
figurations of a microscopic height z = [zi, Z2, 23) which 
defines a two dimensional interface in five dimensions; 
the black and grey loops are contour lines of projec- 
tions of z along different directions. The microscopic 
heights are defined at the centers of the plaquettes of C 
in the following way. The difference in z between two 
neighboring plaquettes depends on the flavor and direc- 
tion of the bond that they share, and is given by one of 
four vectors (also referred to as colors): A = (— 1, 1, 1), 
B = (1,1,-1), C = (-1,-1,-1), or D = (1,-1,1); 
Fig. [l]d. This height mapping is the same as the one for 
the four-coloring model on the square lattice jl3| . 

Third, we coarse grain the microscopic height. We 
view the configurations of z as consisting of flat domains 
over which z is averaged so as to obtain a coarse grained 
height h |13|. The partition function of the loop model 
which incorporates only the large scale fluctuations of 
h can be written as a functional integral Jl?hexp(— S), 
where S is the effective action. S can be broken up into 
three terms, 



S : 

S h = 



S e + Sb 
K 
~2 



S w 

d 2 x [(c\h) 2 



(d 2 h) 



^ (E • h)R 



d 2 x u>(h) 



(3) 



each having a simple geometrical interpretation: 
i) S e describes the elastic fluctuations of the interface. 
Its form, namely the fact that the elasticity is given by a 
single stiffness constant K, is fixed by the symmetries of 
the loop model: four-fold rotations of the lattice, trans- 
lations in height space, and cyclic permutations of the 
colors HQ. 
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ii) Sb is the above mentioned boundary term that cor- 
rects the weights of loops that wind around a cylinder 
0). R is the scalar curvature, which for a particular 
choice of coordinates on the cylinder is the difference 
of two delta functions centered at the two ends. The 
term Sb therefore has the effect of placing vertex opera- 
tors exp(±iEo • h) at the boundaries of the cylinder. In 
the Coulomb gas representation vertex operators are as- 
sociated with electric charges, and E is referred to as 
the background charge j^]. In the chosen normalization 
for the color-vectors E = (— ne , 0, 0). 

ensures the correct weight of loops in the bulk. 
Namely, as mentioned earlier, the loops of the FPL 2 
model are contours of h, and if the only bulk term in S 
were S e then the two loop orientations would be equiv- 
alent - S e equally weights a step up and a step down in 
height. This is inconsistent with and an additional 

bulk term is necessary. 

It is important to note that a similar construction, 
leading to S e and Sb, has been used previously by many 
authors, and it goes under the name: Coulomb gas ap- 
proach to critical phenomena 0,|l5|]. The important dif- 
ference here is the inclusion of the loop- weight term S w . 
The presence of this term in the effective action fixes the 
value of K, which in the traditional Coulomb gas ap- 
proach would remain unknown, only to be determined 
by an exact value of some critical exponent, typically de- 
rived from a Bethe ansatz solution of the model. For the 
cases where Bethe ansatz solutions do exist (e.g., Potts 
models, 0(n) models) we have checked that the geomet- 
rical approach developed here reproduces t he sam e value 
of the stiffness 



10 



Microscopically, the operator w(h(x)) in S w gener- 
ates the vertex weights A(x) in Fig. Ffld: A = exp(— w). 



[ "renormalized coupling" 
hfx 

W 

This operator can be written as u>(h(x)) = Eo • Q(x) 
f[4|l , where Q(x) is the cross-staggered operator defined 
in Ref. it is a vector valued function of the colors 
around the site x e C. Q(x) is periodic in height space 
and it can be expanded in a Fourier series of vertex oper- 
ators exp(iG-h(x)); the vectors G lie in the lattice which 
is reciprocal to the lattice of height periods, and they are 
the electric charges in the Coulomb gas representation. 
The magnetic charges on the other hand form the lattice 
of height periods, and they are associated with vortex 
configurations of the height with topological charge b 

In the long- wavelength limit we only keep the most rel- 
evant vertex operator appearing in the Fourier expansion 
of Q(x); this is the operator with the smallest scaling di- 
mension. The scaling dimension of a general operator 
with an electro- magnetic charge [G, b] is given by S: 



x[G,b] 



1 



G (G - 2E ) 



K 



(4) 
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Therefore, the loop weight term in Eq. (|3j) can be written 
as: 



S w = Ce J exp(iG Q • h(x)) d 2 x 



where C is a constant independent of n, and Gq = 
(— 7r,0, 7r) minimizes x[G,0]. 

For n < 2 the effective action 5, with S w given by 
Eq. ([5]), describes a Liouville field theory with imagi- 
nary couplings; w(h) is the so-called Liouville poten- 
tial. In order for S to describe a conformal field the- 
ory it is necessary that w(h) is a marginal operator, i.e., 
x[G Q ,0] = 2 §. Moreover, if S is a CFT then w(h) is 
exactly marginal, and the coupling Ceo does not flow un- 
der renormalization. The exactly marginal operator w(h) 
gives rise to a line of fixed points, which is the aforemen- 
tioned critical phase of the loop model. For n > 2 the 
parameter eo in Eq. (|^) becomes pure imaginary. From 
Eq. (^), and the assumption that K increases with n, it 
follows that in this case the scaling dimension of w(h) 
is necessarily less then two - it is relevant. Therefore, 
under renormalization n — > 00, and the loop model is in 
the ordered phase. 

In the Coulomb gas representation of two-dimensional 
critical models the Liouville potential defines the so- 
called screening charges, which were originally intro- 
duced on formal grounds - to ensure the non-vanishing 
of the four-point correlation functions || . In the case of 
loop models we have uncovered a geometrical interpre- 
tation of the Liouville potential: it enforces the correct 
weighting of oriented loops by making the two orienta- 
tions incquivalent. 

The assumption that 10(h) is exactly marginal - which 
we refer to as the conformal ansatz - has a simple ge- 
ometrical interpretation. Namely, the loop weight (n = 
2cos(-7reo)) is thermodynamically conjugate to the num- 
ber of loops (N), and its non-renormalizability implies 
that the number of large loops, which are responsible for 
the long- wavelength fluctuations of the height, does not 
change under renormalization. In this form the conformal 
ansatz is closely analogous to the hypothesis often made 
for critical percolation, that the number of spanning clus- 
ters is of order one, which in turn leads to hyperscaling 

c. Exact results The conformal ansatz implies 
i[Gq,0] = 2, and from Eq. (^J) we calculate the exact 
value of the stiffness 



K = -(l- e ) 



(6) 



for further convenience we define the parameter g = 
AK/tt. Using Eqs. (g) and (||) for the coupling constants 
of the effective field theory (Eq. (|3|)), we can calculate the 
conformal charge and the critical exponents of the loop 
model along the critical line n < 2. 
The conformal charge is pf: 



c = 3 - 12 
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= 3-12 



(1-.9) 2 



(7) 



(5) 



For n = 2 and n = 1, g = 1 and g — 2/3, and conse- 
quently c = 3 and c = 1; c = 3 is the conformal charge 



3 



of the four-coloring model which can be mapped to the 
n = 2 FPL 2 model (b|], while c = 1 is the known con- 
formal charge of the equal- weighted six- vertex model , 
which is equivalent to the n = 1 FPL 2 model . 

In the limit of Hamiltonian walks (n — > 0, g — 1/2) we 
find c = — 3. This differs from the conformal charge for a 
single Hamiltonian walk on the honeycomb (c = — 1 |]]) 
or the square (c = —1.00(1) |ll|| ) lattice, and we conclude 
that the FPL 2 model defines a new universality class of 
self-avoiding random walks. This is also confirmed by the 
critical exponents, which we calculate next. 

The most commonly studied operators in the con- 
text of loop models are the so called string ("water- 
melon") operators whose two-point correlation function 
gives the probability of having m loop segments propa- 
gating between two points on the lattice [Q~7| ; here we fo- 
cus on string operators associated with black loops only. 
In the Coulomb gas representation of the FPL 2 model 
these operators are associated with mixed electric and 
magnetic charges: [(— e 7r/2, 0, eo7r/2), (— 2k, 0, — 2k)] is 
the charge for the m = 2k string operator, while 
[(-e 7r,0,0),(-2fc,-2,2fc-2)]form = 2fe-l Q. These 
charges can be determined microscopically by mapping 
the string configurations to vortex configurations of the 
height; an analogous analysis for the FPL model on the 
honeycomb lattice can be found in Ref. jio) . From Eq. 
we can calculate the scaling dimensions of the string op- 
erators, 



9,2 (1 - 9? 
X2k =2 k 2g~ 

«^ = §(*»- * + 

2 9 



(8) 



The dimensions associated with an even number of 
strings are the same as those found for the FPL model on 
the honeycomb lattice, while the odd ones are new. In 
particular, in the Hamiltonian limit (n — ► 0, g = 1/2) we 
find xi = —1/4. This, via standard scaling relations M, 
leads to a prediction for the exponent 7 = 1 — x\ = 5/4, 
which describes the scaling of the number of walks with 
the number of steps (bonds traversed). This value of 7 
is different from the ones found for a single walk on the 
honeycomb (7 = 1 ||) and the square lattice (7 s» 1.0444 
jllj ) , confirming once again that the FPL 2 model defines 
a new universality class of Hamiltonian walks. 

For n = 1 the FPL 2 model coincides with the FPL 
model of Batchelor et al. [jll). For this value of n, 
g = 2/3, and x% — 1/6 and xi = 1/4, while the nu- 
merical transfer matrix results quoted in Ref. p"l| ] are 
xi = 0.1667(1) and x 2 = 0.2500(1). Furthermore Batch- 
elor et al. found numerical evidence for the relation 
X3 = x± + g, which we find to be satisfied exactly. 

In conclusion, we have shown that the conformal field 
theory associated with the scaling limit of a lattice loop- 
model can be constructed by mapping it to a random 
surface. The loops appear as contour lines of a surface 
whose long wavelength fluctuations are described by a 



Liouville field theory. The remarkable feature of loop 
models is that the couplings in the field theory are com- 
pletely fixed by geometry thus leading to exact results 
for their critical properties. 

The conformal ansatz, which is equivalent to the state- 
ment that the number of large loops does not change with 
length scale, is central to the whole construction, and it 
leads to a physical interpretation of the Dotsenko-Fateev 
screening charges [[| . The validity of this ansatz has been 
confirmed for loop models for which an exact solution ex- 
ists (e.g., Potts models, 0(n) models) Q,|o|. In order to 
show this explicitly a renormalization group which takes 
into account the extended nature of loops will have to be 
developed. This remains an interesting open question. 
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